Chameleon dark energy models with characteristic signatures 
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In chameleon dark energy models, local gravity constraints tend to rule out parameters in which 
observable cosmological signatures can be found. We study viable chameleon potentials consistent 
with a number of recent observational and experimental bounds. A novel chameleon field potential, 
motivated by f(R) gravity, is constructed where observable cosmological signatures are present both 
at the background evolution and in the growth-rate of the perturbations. We study the evolution 
of matter density perturbations on low redshifts for this potential and show that the growth index 
today 70 can have significant dispersion on scales relevant for large scale structures. The values of 70 
can be even smaller than 0.2 with large variations of 7 on very low redshifts for the model parameters 
constrained by local gravity tests. This gives a possibility to clearly distinguish these chameleon 
models from the A-Cold-Dark-Matter (ACDM) model in future high-precision observations. 

PACS numbers: 04.50.Kd, 95.36.+X 



I. INTRODUCTION 

The accelerated expansion of the Universe today is a 
very important challenge faced by cosmologists For 
an isotropic comoving perfect fluid, a substantially nega- 
tive pressure is required to give rise to the cosmic accel- 
eration. One of the simplest candidates for dark energy 
is the cosmological constant with an equation of state 
wde = — 1) but we generally encounter a problem to ex- 
plain its tiny energy density consistent with observations 

There are alternative models of dark energy to the 
cosmological constant scenario. One of such models is 
quintessence based on a minimally coupled scalar field 
with a self- interacting potential Q. In order to realize 
the cosmic acceleration today, the mass of quintessence 
is required to be very small (m^ « 10 -33 GeV). From a 
viewpoint of particle physics, such a light scalar field may 
mediate a long range force with standard model particles 
For example, the string dilaton can lead to the vio- 
lation of equivalence principle through the coupling with 
baryons In such cases we need to find some mecha- 
nism to suppress the fifth force for the consistency with 
local gravity experiments. 

There are several different ways to screen the field in- 
teraction with baryons. One is the so-called run away 
dilaton scenario [6j in which the field coupling F((f>) 
with the Ricci scalar R is assumed to approach a con- 
stant value as the dilaton </> grows in time (e.g., F(<j)) = 
C\ + C-2e^^ as (j> — > 00). Another way is to consider a 
field potential having a large mass in the region of high 
density where local gravity experiments are carried out. 
In this case the field does not propagate freely in the local 
region, while on cosmological scales the field mass can be 
light enough to be responsible for dark energy. 

The latter scenario is called the chameleon mechanism 



in which a density-dependent matter coupling with the 
field can allow the possibility to suppress an effective cou- 
pling between matter and the field outside a spherically 
symmetric body 0, H| ■ The chameleon mechanism can 
be applied to some scalar-tensor theories such as f{R) 
gravity H, [Io[ and Brans-Dicke theory [TO]. In f(R) 
gravity, for example, there have been a number of viable 
dark energy models [l2| that can satisfy both cosmolog- 
ical and local gravity constraints. For such models the 
potential of an effective scalar degree of freedom (called 
"scalaron" [l| ) in the Einstein frame is designed to have 
a large mass in the region of high density. Even with 
a strong coupling between the scalaron and the baryons 
(Q = — 1/V6), the chameleon mechanism allows the f(R) 
models to be consistent with local gravity constraints. 

The chameleon models are a kind of coupled 
quintessence models [l4[ defined in the Einstein frame 
0, Q. While the gravitational action is described by 
the usual Einstein-Hilbert action, non-rclativistic mat- 
ter components are coupled to the Einstein frame metric 
multiplied by some conformal factor which depends on 
a scalar (chameleon) field. This is how the gravitational 
force felt by matter is modified. While there have been 
many studies for experimental and observational aspects 
of the chameleon models [H|-[38j, it is not clear which 
chameleon potentials are viable if the same field is to be 
responsible for dark energy. 

In this paper we identify a number of chameleon po- 
tentials that can be consistent with dark energy as well as 
local gravity experiments. We then constrain the viable 
model parameter space by using the recent experimental 
and observational bounds- such as the 2006 Eot-Wash 
experiment (39j . the Lunar Laser Ranging experiment 
[40j and the WMAP constraint on the time- variation of 
particle masses [41[ . This can actually rule out some of 
the chameleon potentials with natural model parameters 
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for the matter coupling Q of the order of unity. 

In order to distinguish the viable chameleon dark en- 
ergy models from the ACDM model, it is crucial to study 
both the modifications in the evolution of the background 
cosmology and the modified evolution of the cosmological 
density perturbations. For the former, we shall consider 
the evolution of the so-called statefindcrs introduced in 
Refs. [42J, |43| and show that these parameters can exhibit 
a peculiar behavior different from those in the ACDM 
model. 

On the other hand, the growth "index" 7 of matter per- 
turbations d m defined through dln5 m /dlna = (f2^) 7 , 
where a is a scale factor and Q*^ is the density param- 
eter of non-relativistic matter, is an important quantity 
that allows to discriminate between different dark energy 
models and interest in this quantity was revived in the 
context of dark energy models [3 E|| . 

Its main importance for the study of dark energy mod- 
els stems from the fact that for ACDM the quantity 7 is 
known to be nearly constant with respect to the red- 
shift z, i.e. 7acdm = 7o — 0.02,2 to exquisite accuracy 
HI], with 70 = j(z = 0) « 0.555 [3. As emphasized 
in Ref. [46j . large variations of 7 on low redshifts could 
signal that we are dealing with a dark energy model out- 
side General Relativity. This was indeed found for some 
scalar-tensor dark energy models |47| and f(R) models 
[48M50l ] . Such large variations can also occur in models 
where dark energy interacts with matter [5lL [52j . This is 
exactly the case in chameleon models, which we investi- 
gate in this paper, because of the direct coupling between 
the chameleon field <j> and all dust-like matter. This di- 
rect coupling is however not confined to the dark sector 
as in standard coupled quintessence. 

An additional important point is the possible appear- 
ance of a scale-dependence or dispersion in 7. Hence 
the behavior of 7 on low redshifts can be both time- 
dependent and scale-dependent j53l - l55| . This dispersion 
can also be present in the models investigated here. Using 
the observations of large scale structure and weak lensing 
surveys, one can hope to detect such peculiar behaviors 
of 7 (see e.g., Ref. [56|). If this is the case this would sig- 
nal that the gravitational law ma y b e modified on scales 
relevant to large scale structures (45l IBH l55l [571462^ . 

In this paper we study the evolution of 7 as well as its 
dispersion, its dependence on the wavenumbers of pertur- 
bations. We shall show that some of the chameleon mod- 
els investigated here can be clearly distinguished from 
ACDM through the behavior of 7 exhibiting both large 
variations and significant dispersion, with the possibility 
to obtain small values of 7 today as low as 70 < 0.2. 



II. CHAMELEON COSMOLOGY 

A. Background equations 

In this section we review the basic background evolu- 
tion of chameleon cosmology. We consider the chameleon 



theory described by the following action 0, H, [63| 



S = J d A x^g~ 



+S m [H m ;A 2 {(t>) 



(1) 



where g is the determinant of the (Einstein frame) met- 
ric g^, R is the Ricci scalar, G is the bare gravitational 
constant, is a scalar field with a potential V((f>), and 
S m is the matter action with matter fields At low 
redshifts it is sufficient to consider only non-relativistic 
matter (cold dark matter and baryons), but for a gen- 
eral dynamical analysis including high redshifts radiation 
must be included. 

We assume that non-relativistic matter is universally 
coupled to the (Jordan frame) metric A 2 ((f>) c/ M „, the Ein- 
stein frame metric g^ v multiplied by a field-dependent 
(conformal) factor A 2 (<fr). This direct coupling to the 
field 4> is how the gravitational interaction is modified. 
We can generalize this to arbitrary functions A^> (cf>) for 
each matter component pi, but in this work we will take 
the same function A((f>) for all components. We write the 
function A((f>) in the form 



A{4>) 



(2) 



where M p i = 1/ V8ttG is the reduced Planck mass and Q 
describes the strength of the coupling between the field 
</> and non-relativistic matter. In the following we shall 
consider the case in which Q is constant. In fact, the con- 
stant coupling arises for Brans-Dicke theory b y a confor- 
mal transformation to the Einstein frame [§l lllj. Even 
when \Q\ is of the order of unity, it is possible to make 
the effective coupling between the field and matter small 
through the chameleon mechanism. 

Let us consider the scalar field cf> together with non- 
relativistic matter (density p^) and radiation (density 
p r ) in a spatially flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) space-time with a time-dependent scale 
factor a(t) and a metric 



ds 2 = g^ y dx^ dx v 



-dt 2 



\t) dx 2 



(3) 



The corresponding background equations are given by 
3H 2 = ( n + p* m + p r ) /M 2 , , (4) 
+ 3Hj> + V,t = -Qp m /M pU (5) 
p* m + 3Hp* m = Q P ;J/M pl , (6) 

p r + AHpr = , (7) 

where p^ = <fi 2 /2 + V{4>), = dV/d<p, and a dot rep- 
resents a derivative with respect to cosmic time t. The 
quantity p* n is the energy density of non-relativistic mat- 
ter in the Einstein frame and we have kept the star to 
avoid any confusion. Integration of Eq. (|6|) gives the so- 
lution p* m oc a _3 e ( 5^/^ f p 1 . We define the conserved matter 
density: 



-Q<t>/M p i 



(8) 
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which satisfies the standard continuity equation, p m + 
3Hp m = 0. Then the field equation ([5]) can be written 
in the form 



ZH<t> + V cS .4, = , 



(9) 



where V e s is the effective potential defined by 

V oS = V(4>) + e^'ft, . (10) 

We emphasize that it is the Einstein frame which is the 
physical frame. Due to the coupling between the field 
and matter, particle masses do evolve with time in our 
model. 

We consider runaway positive potentials V(<fi) in the 
region <f> > 0, which monotonically decrease and have a 
positive mass squared, i.e. < and > 0. We 
also demand the following conditions 



lim 

0->O 



Yd 
v 



lim 



V, 



0— )-oo V 



0. 



(11) 



The former is required to have a large mass in the region 
of high density, whereas we need the latter condition to 
realize the late-time cosmic acceleration in the region of 
low density. At <fi = the potential approaches either 
oo or a finite positive value Vq. In the limit <p — > oo wc 
have either V — > or V — > Voo, where Voo is a nonzero 
positive constant. If Q > the effective potential V c s((j)) 
has a minimum at the field value cf> m (> 0) satisfying the 
condition V e s^(4> m ) = 0, i.e. 



V^ m ) + Q(p m /M pl )eQ^^ = . 



(12) 



If the potential satisfies the conditions > and 
V,4>4> > in the region </> < 0, there exists a minimum at 
4> = 4>m (< 0) provided that Q < 0. In fact this situation 
arises in the context of f(R) dark energy models [9l Hoj|. 
Since the analysis in the latter is equivalent to that in 
the former, we shall focus on the case Q > and < 
in the following discussion. 



B. Dynamical system 

In order to discuss cosmological dynamics, it is conve- 
nient to introduce the following dimensionless variables 



V 



V6HM pl ' 



X2 



X3 



pi 



(13) 

Equation (|4]) expresses the constraint existing between 
these variables, i.e. 



where 



2 _,_ 2 

X l "T X 2 ! 



1 A\ j-- 



(14) 



(15) 



Taking the time-derivative of Eq. (J4j) and making use of 
Eqs. ([5])- ([7]), it is straightforward to derive the following 
equation 



(3 + 3xj - 3a 



(16) 



where a prime represents a derivative with respect to 
N = In a. A useful quantity is the effective equation of 
state 



Wcff = -1 - = x\ - X% + ^ 



(17) 



We also introduce the field equation of state as 
2 /2 - V{4>) _ x\ - xl 



2 /2 + V{4>) xl + : 



(18) 



Using Eqs. ([!])- 0, we obtain the following equations 



V6. 



H' 



x 1 = —3xi H — tt^x 2 — Xi — 



H 



V6 
2 
V6 



Q (l x ^ x 2 £3) , 
H' 



x' 2 = — — Xx 1 x 2 - x 2 — 
1 rl 

H' 

X 3 = -2x 3 - X 3^j > 

X' = -V&X 2 (T-l)xi, 



where 



A = — 



p\V,<f, 



V 



r = 



vv. 



(19) 
(20) 

(21) 
(22) 

(23) 



From the conditions (fTT|) it follows that the quantity A 
decreases from 00 to as <f> grows from to 00. Since 
x\ > in Eq. (|22|) . the condition A' < translates into 



vl 



(24) 



Chameleon potentials shallower than the exponential po- 
tential (r = 1) can satisfy this condition. 

Once the field settles down at the minimum of the ef- 
fective potential (Tl7J|) , we have 



x x ~ 0, 



X2 



A 



(25) 



which gives ~ — 1 from Eq. (|18p . As the matter den- 
sity p* m decreases, the field evolves slowly along the in- 
stantaneous minima characterized by ([25]). Wc require 
that A Q = O(l) during radiation and deep matter 
eras for consistency with local gravity constraints in the 
region of high density. For the dynamical system (|19|) - 
(|2"2"j) there is another fixed point called the "(^-matter- 
dominated era (0MDE)" (Tj where fl^ = w cS = 2Q 2 /3. 
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However, since we are considering the case in which Q is 
of the order of unity, the effective equation of state w e g is 
too large to be compatible with observations. Only when 
Q < 0(0.1) the 0MDE can be responsible for the matter 
era [H|. 

When the chameleon is slow-rolling along the mini- 
mum, we obtain the following relation from Eqs. (|15p 
and ([25]): 



Q n 4 



(26) 



While A Q during the radiation and matter eras, A 
becomes the same order as Q around the present epoch. 
The field potential is the dominant contribution on the 
r.h.s. of Eq. (J4| today, so that 



V{<fo) - 3H£M^ ~ Pc , 



(27) 



where the subscript "0" represents present values and 
p c ~ 10~ 29 g/cm 3 is the critical density today. 

III. CHAMELEON MECHANISM 

In this section we review the chameleon mechanism as 
a way to escape local gravity constraints. In addition 
to the cosmological constraints discussed in the previ- 
ous section, this will enable us to restrict the forms of 
chameleon potentials. 

Let us consider a spherically symmetric space-time in 
the weak gravitational background with the neglect of the 
backreaction of metric perturbations. As in the previous 
section we consider the case in which couplings Qi are the 
same for each matter component (Qi = Q), i.e., in which 
the function A(<fi) is given by Eq. (JSJ). Varying the action 
(jTJ) with respect to (f> in the Minkowski background, we 
obtain the field equation 



d 2 4> 

dr 2 



2d<j) 
r dr 



dV cS 
d<f> 



(28) 



where r is the distance from the center of symmetry and 
V e g is defined in Eq. (fTUj) . 

Assuming that a spherically symmetric object (radius 
r c and mass M c ) has a constant density p m = pa with 
a homogeneous density p m = pg outside the body, the 
effective potential has two minima at <fr = 4>a and <j> = <j>g 
satisfying the conditions 



v 4 

Va 



Q(p A /M p} )e^ M * =0, (29) 
Q(p B /M p i)e Q + B ' M * = 0. (30) 



Since Q4*a/M p \ <C 1 and Q4>b/M p \ <C 1 for viable field 
potentials in the regions of high density, the conserved 
matter density p m is practically indistinguishable from 
the matter density in the Einstein frame. 

The field profile inside and outside the body can be 
found analytically. Originally this was derived in Refs. 0, 



|8j under the assumption that the field is frozen around 
4> = 4>a in the region < r < r±, where r\ (< r c ) is the 
distance at which the field begins to evolve. It is possible, 
even without this assumption, to derive analytic solutions 
by considering boundary conditions at the center of the 
body [26[. 

We consider the case in which the mass squared 
m 2 B = d d ^ ff ((f) b) outside the body satisfies the condi- 
tion rriBT c <C 1, so that the m^-dependent terms can be 
negligible when we match solutions at r — r c . The re- 
sulting field profile outside the body (r > r c ) is given by 

M 



4,(r) = (f>i 



Qcs M c 



4irM, 



pi 



(31) 



where the effective coupling Qeff between the field and 
matter is 



Q 



r° r c (m A r c ) 
mAT\(e mAri + e~ 



(32) 



The mass rriA is defined by m 2 A = d d ^ ff (<Pa)- 
The distance r\ is determined by the condition 
171 A [ < ^( r i) — 4>a] = Qpa, which translates into 



*B- <t>A + Qp A (rl-K)l(2M vi ) 
_ 6QM p i<Z> c m A ri(e mAri + e~ mAri ) 
' (m A r c ) 2 e mATl - e~ m * r i 



(33) 



where <& c = M c /(%nr c ) = pat 2 /(6M 2 ^) is the gravita- 
tional potential at the surface of the body. 

The fifth force exerting on a test particle of a unit mass 
and a coupling Q is given by F = —QW<j)/Mp\. Using 
Eq. ([BT]) . the amplitude of the fifth force in the region 
r > r c is 



F 



2\QQ cS \^ 



(34) 



As long as | Q c g | <C 1 , it is possible to make the fifth force 
suppressed relative to the gravitational force GM c /r 2 . 



From Eq. (|32|) the effective coupling Q c g can be made 
much smaller than Q provided that the conditions Ar c = 
r c — ?*i <C r c and vriAfc 3> 1 are satisfied. Hence we 
require that the body has a thin-shell and that the field 
is heavy inside the body for the chameleon mechanism to 
work. 

When the body has a thin-shell (Ar c -c r c ), one can 
expand Eq. (|33[) in terms of the small parameters Ar c / r c 
and l/(m^r c ). This leads to 



6QA/ P i$ c 



Arc 



1 



m A r c 



(35) 



where eth is called the thin-shell parameter. As long as 
tuaTc S> (A^c/rc)" 1 , this recovers the relation e t h — 
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Ar c /r c 0, Q . The effective coupling (l32j) is approxi- 
mately given by 



Q e g ~ 3Qe t h ■ 



(36) 



If eth is much smaller than 1 then one has Q c s <C Q, 
so that the models can be consistent with local gravity 
constraints. 

As an example, let us consider the experimental bound 
that comes from the solar system tests of the equiva- 
lence principle, namely the Lunar Laser Ranging (LLR) 
experiment, using the free-fall acceleration of the Moon 
(aMoon) and the Earth (a©) toward the Sun (mass M Q ) 
[a H3, [26| ■ The experimental bound on the difference of 
two accelerations is given by 



2|omc 



a© | 



(«Mq 



«©) 



< 10 



-13 



(37) 



Under the conditions that the Earth, the Sun, and the 
Moon have thin-shells, the field profiles outside the bod- 
ies are given as in Eq. (|3ip with the replacement of corre- 
sponding quantities. The acceleration induced by a fifth 
force with the field profile (f>(r) and the effective coupling 
Q off is a fifth = \Q eS V<f>(r)/M v i\. Using the thin-shell pa- 
rameter £tx^ for the Earth, the accelerations and 
aMoon are 



GAL 







r 

GM, 



l + 18Q 2 4, e |p 







1 



18Q 2 e. 2 



<I>- 



th,ffi 



O^Moon 



(38) 
(39) 



where $© ~ 2.1 x 1CT 6 , $© ~ 7.0 x 10~ 10 , and $ M oon ^ 
3.1 x 10~ n are the gravitational potentials of Sun, Earth 
and Moon, respectively. Then the condition (|3T|) reads 



Cth,0 < 



8.8 x 10 



-7 



Q 



(40) 



Using the value $ ffi ~ 7.0 x 10~ 10 , the bound (J3DJ trans- 
lates into 



^,© < lO-^M, 



(41) 



where we used the condition 4>b,® ^ 4>A.@- For the Earth 
one has pa — 5g/cm 3 (mean density of the Earth) 3> 
Pb — 10~ 24 g/cm 3 (dark matter/baryon density in our 
galaxy), so that the condition <pB,@ 3> 'Pa.® is satisfied. 

In Sec. lIVI we constrain viable chameleon potentials by 
employing the condition (f4"T| together with the cosmo- 
logical condition we discussed in Sec. [TTJ In Sec. [V] we 
restrict the allowed model parameter space further by 
using a number of recent local gravity and observational 
constraints. 



cosmological constraints. Let us consider the potential 

V(4>) = A/ 4 /(</>) , (42) 

where M is a mass scale and f(4>) is a dimensionlcss 
function in terms of <f>. 

The local gravity constraint coming from the LLR ex- 
periment is given by Eq. (|41[) . where <Pb,® is determined 
by solving 



{mm^b^I-Qpb/m 4 . 



(43) 



Here we take ps — 10~ 24 g/cm 3 for the homogeneous 
density outside the Earth. Once the form of /(</>) is spec- 
ified, the constraint on the model parameter, e.g., M, can 
be derived. 

From the cosmological constraint (f26|) . \/Q is of the 
order of 1 today. Then it follows that 



|Mpi/, (^,)| ~ g/(0 o ) * QPc/M A , 



where we used Eq. (|27|) . 

We also require the condition (j2~4"|) , i.e. 

y, _ f f,4>4> , 
1 f2 ^ 1 ' 

J .A 



(44) 



(45) 



for all positive values of <fi. We shall proceed to find viable 
potentials satisfying the conditions (T4"Tj) . (|4"3"1) . (f4"4"l) , and 
(|45|) . From Eqs. (g3J) and rtHJ) we obtain 



J, 



PB 



JAM 



10 5 



(46) 



Let us consider the inverse power-law potential V((j)) 
M 4+n cj)- n (n > 0), i.e. 



m = (M/<t>y 



(47) 



Since T = (n + l)/n > 1, the condition (|45|) is automati- 
cally satisfied. The cosmological constraint (j4~4")) gives 

<fo ~ nM pl /Q . (48) 

From Eq. (|46|) we find the relation between <f>o and 4>b,® '■ 

fa ~ 10 5 /(" +1 V B , e . (49) 

Using the LLR bound (|4T|) . it follows that 

(50) 



5(3n + 2) 
K < 10 »+» Mpl, 



This is incompatible with the cosmological constraint 
(|4"8"f for n > 1 and Q = 0(1)- Hence the inverse power- 
law potential is not viable. 



IV. VIABLE CHAMELEON POTENTIALS 



A. Inverse power-law potential + constant 



We now discuss the forms of viable field potentials that 
can be in principle consistent with both local gravity and 



The reason why the inverse power-law potential does 
not work is that the field value today required for cosmic 
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acceleration is of the order of M p i, while the local gravity 
constraint demands a much smaller value. This problem 
can be circumvented by taking into account a constant 
term to the inverse power-law potential. Let us then 
consider the potential V{<j>) = M 4 [1 + fj,(M/(ft) n ] (n > 
0), i.e. 



f{4>) =i + ^M/4>y 



(51) 



where u is a positive constant. The rescaling of the mass 
term M always allows to normalize the constant to be 
unity in Eq. (|5ip . For this potential the quantity T reads 



1 



1 



1 



(52) 



which satisfies the condition T > 1. In the region 
/i(M/(f)) n 3> 1 we have that V ~ (n+l)/n, which recovers 
the case of the inverse power-law potential. Meanwhile, 
in the region /j,(M/<fi) n <§C 1, one has r> 1. The latter 
property comes from the fact that the potential becomes 
shallower as the field <f> increases. This modification of 
the potential allows a possibility that the model can be 
consistent with both cosmological and local gravity con- 
straints. 

The addition of a constant term to the inverse power- 
law potential does not affect the condition (|46| . which 
means that the resulting bounds (|49[) and ([5T)|) are not 
subject to change. On the other hand, the cosmologi- 
cal constraint (|48|) is modified. Let us consider the case 
where the condition fi(M/<fio) n 1 is satisfied today, i.e. 
f ((t) a ) ~ 1. From Eq. (JUJ) it follows that 



M 



P, 



1/4 



10" 12 GeV. 



(53) 



and 



l/(n+l) 



nu 



M, 



M p i ~ (10 



-30n 



(54) 

Hence the field value </>o today can be much smaller than 
the Planck mass, unlike the inverse power-law potential. 
From Eqs. (f50|) and (|54|) we get the constraint 



< 10 



15n-10 



In. 



(55) 



If n = 1, for example, one has a < 10 5 . For 
larger n the bound on fi becomes even weaker. Wc 
note that the condition n(M/<f>) n < 1 is satisfied for 
0/Mpi > 10" 10 /"- 15 /" 17 ™- This shows that even the 
field value such as 4>b.@ = 10~ 15 Af p i satisfies the con- 
dition fj,(M/(f) BiS ) n < 1. Thus the term a{M/(j)) n is 
smaller than 1 for the field values we are interested in 

{4>b,® < </> < M- 

A large range of experimental bounds for this model 
has been derived in the literature, see Refs. [1, [l?], Ho|, H3 ■ 
For Q = n = 1 . it was found in Ref . [2(| that the model is 
ruled out by the Eot-Wash experiment unless fi < 10~ 5 . 
This applies for general n: to obtain a viable model for 



Q of the order of unity one must impose a fine-tuning 
n 3> 1 or n <C 1. 

The potentials, which have only one mass scale equiv- 
alent to the dark energy scale, are usually strongly con- 
strained by the Eot-Wash experiment. We shall look into 
this issue in more details in Sec. fVl 



B. Construction of viable chameleon potentials 
relevant to dark energy 

The discussion given above shows that a function 
/(</>) that monotonically decreases without a constant 
term is difficult to satisfy both cosmological and lo- 
cal gravity constraints. This is associated with the 
fact that for any power-law form of f((j>) the condition 
\M p if y( j,((j)o)/f{(j)o) \ — 1 leads to the overall scaling of the 
function /(0o) itself, giving </>o of the order of M p \. The 
dominance of a constant term in f(4>o) changes this sit- 
uation, which allows a much smaller value of 0o relative 
to M p i. 

Another example similar to V(4>) = M 4 [l + /j,(M/<j>) n ] 
is the potential [13 

F(^)=M 4 exp[M(M/^r], (56) 

where fi > and n > 0. For this model the quantity 

r_i + 2±iI(*Y\ (57, 

is larger than 1. In the asymptotic regimes characterized 
by u(M/<j)) n > 1 and u(M/(/)) n < 1 we have r ~ 1 and 
T > 1, respectively. When u(M/4>) n < 1 the function 
/((/>) = exp[u(M / <f)) n ] can be approximated as /((/>) ~ 
1 + fj,(M/(f)) n , which corresponds to Eq. (f5T|) . In this case 
the constraints on the model parameters are the same as 
those given in Eqs. ([53"!) - ([55)) . 

There is another class of potentials that behaves as 
V{(f) ~ A/ 4 [l - n(<f>/M) n ] (n > 0,0 < n < 1) in the re- 
gion n((j)/M) n -C 1. In fact this asymptotic form corre- 
sponds to the potential that appears in f(R) dark energy 
models. While the potential is finite at <\> = the deriva- 
tive \V$\ diverges as (j) —¥ for < n < 1, so that the 
first of the condition (fTTj) is satisfied. In order to keep 
the potential positive we need some modification of V in 
the region ^{(j)/M) n > 1. 

In scalar-tensor theory it was shown in Ref. [ll| that 
the Jordan frame potential of the form U(<p) = i\/ 4 [l — 
fi(l - e - 2< 3*/ M p')™] (0 < n < 1, < n < 1) can satisfy 
both cosmological and local gravity constraints. In this 
case the potential V(4>) in the Einstein frame is given 
by V((f>) = e 4< 5^/ Mpl L r (0), which possesses a de Sitter 
minimum due to the presence of the conformal factor. 
Cosmologically the solutions finally approach the de Sit- 
ter fixed point, so that the late-time cosmic acceleration 
can be realized. 

Now we would like to consider a runaway positive po- 
tential in the Einstein frame. One example is 

-0/M p i\Tll 



V((j>) = M 4 [l - n(l - e 



(58) 
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FIG. 1: The potential ([58]) versus the field tf> for fi = 0.7 and 
n = 0.7. The potential has a finite value V = Af 4 at = 0, 
but its derivatives diverge (|V^| — > oo and — > oo) as 



where < /i < 1 and < n < 1. 
haves as V(cj>) ~ Af 4 [l - /i(0/M p i) 7 
approaches V(0) — >• Af 4 (l — /i) in the limit 
Fig. [J). For the potential flSSJ) we obtain 



This potential bc- 
for 4> < A/pi and 



» M p i (see 



r - 1 = 



1 — nx — fJ.(l — x) r 
n/ix(l — x) n 



(59) 



One can easily show that the r.h.s. is positive under the 
conditions 0</i<l, 0<n<l, so that T > 1. In both 
limits 0—^0 and (j> —> oo one has T — > +oo. Since T has 
a minimum at a finite field value, the condition T ^> 1 is 
not necessarily satisfied today (unlike the potential ([56])). 

Unless is very close to 1 the potential energy today 
is roughly of the order of M 4 , i.e. f(<po) ~ 1- From 
Eq. ((331) it; tnen follows that 



- x Q ) n 1 x ~ Q , 



(60) 



where x = e - ^ ^ 1 . If O < M p i, we have that 
<j>o/M p i ~ (n/i) 1 -^ 1- "). From Eq. (g3J) we obtain 



npt(l-.T B ) ^ 10 5 Q 



(61) 



where is = e ^ b / m p 1 . Under the condition 4>b <■ M p \ 
we have B /A/ p i ~ (10- 5 n///Q) 1/(1_n) from Eq. f6Tj). 
Then the LLR bound (j41|) corresponds to 



(62) 



10 10-15« < 



When /i = 0.5 and = 0.05 with Q = 1, the constraint 
((6"2"j) gives n > 0.63 and n > 0.56 respectively. 



C. Statefinder analysis 

The statefinder diagnostics introduced in Refs. (42l.l43l| 
can be a useful tool to distinguish dark energy models 
from the ACDM model. The statefinder parameters are 
defined by 



r-1 



iH 3 



3(9-1/2)' 



(63) 



where q = —a/(aH 2 ) is the deceleration parameter. 
Defining h = H 2 . it follows that 
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_ h!_ 



1+ 2h 



3ti 
~2h 



"oTT ' (64) 



where a prime represents a derivative with respect to 
N = In a. 

In the radiation dominated epoch we have h oc e~ 4Ar , 
which gives (r, s) ~ (3, 4/3). During the matter era r ap- 
proaches 1, whereas s blows up from positive to negative 
because of the divergence of the denominator in s (i.e. 
q = 1/2). For the chameleon potentials ([5"6"|) and l|58p 
the solutions finally approach the de Sitter fixed point 
characterized by (r, s) = (1,0). Around the de Sitter 
point the solutions evolve along the instantaneous min- 
ima characterized by (x\, X2, X3) = (X/s/6, \J\ — A 2 /6, 0) 
with h oc V. Using Eq. (f22|) as well, one has h'/h ~ — A 2 
and h" //i ~ (2T— 1)A 4 . Then the statefinder diagnostics 
around the de Sitter point can be estimated as 



1 



2 J 2 



(2F - 1)A 4 - 3A 2 
3(3- A 2 ) 



(65) 
(66) 



Since A finally approaches 0, Eq. ([2"2"| implies that 
TA 3 — > asymptotically. For the potential in which 
r ^> 1 holds today it can happen that TA 4 ^> A 2 , which 
gives r ~ 1 + TA 4 > 1 and s ~ -2rA 4 /[3(3 - A 2 )] < 
around the present epoch. In the upper panel of Fig. [5] 
we plot the evolution of the variables r and s for the po- 
tential (f56|) in the redshift regime —1 < z = 1/a—l < 10. 
The statefinders evolve toward the de Sitter point char- 
acterized by (r, s) = (1,0) from the regime r > 1 and 
s < 0. This behavior is different from quintessence with 
the power- law potential V{(j>) = M i+n cj}~ n ( n > 0) in 
which the statefinders are confined in the region r < 1 
and s>0@|. 

The potential (|58p allows the possibility that T is not 
much larger than 1 even at the present epoch. In the 
regime A 2 <C 1 we then have r ~ 1 — 3A 2 /2 < 1 and 
s ~ A 2 /3 > 0. In fact we have numerically confirmed 
that the solutions enter this regime by today (see the 
lower panel of Fig. [2]). Finally they approach the de Sitter 
point from the regime r < 1 and s > 0. Hence one 
can distinguish between chameleon potentials from the 
evolution of statefinders. 
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FIG. 2: (Top): Evolution of the statefinders r and s for the 
potential (|56|) with n = 1, Q = 1 and fi — 1. The solutions 
approach the de Sitter point at (r, s) = (1, 0) from the region 
r > 1 and s < 0. (Bottom): Evolution of statefinders for the 
potential ()58[) with n = 0.7, Q — 1 and = 0.7. In this case 
the solutions approach the de Sitter point from the region 
r < 1 and s > 0. 
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LOCAL GRAVITY CONSTRAINTS ON 
CHAMELEON POTENTIALS 



In this section we discuss a number of local gravity 
constraints on the chameleon potentials ([55)1 and (|58p in 
details. Together with the LLR bound (|41 [) we use the 
constraint coming from 2006 Eot-Wash experiments [39| 
as well as the WMAP bound on the variation of the field- 
dependent mass. 



A. The WMAP constraint on the variation of the 
particle mass 

Due to the conformal coupling of the field <j> to matter, 
any particle will acquire a 0-dependent mass: 



(0) = moe^/^ 1 



(67) 



where m is a constant. 

The WMAP data constrain any variation in m (</>), be- 
tween now and the e poc h of recombination to be < 5% 
at 2cr (< 23% at 4<r) [4l|. We then require that 



Am(<£) 




m 





3(4>O-0rcc) 



< 0.05, 



(68) 



where rcc is the field value at the recombination epoch. 
If we assume that the chameleon follows the minimum 
since recombination then 4>o 3> r oc, and the field in the 
cosmological background today must satisfy 



Qh/M p i < 0.05 . 



(69) 



This provides a constraint on the coupling Q and the 
model parameters of chameleon potentials. 

Note that the WMAP constraint is not a local gravity 
constraint. Nevertheless, it provides strong constraints 
on the potential ([58]) with natural parameters and is 
therefore considered in this section. 



B. Constraints from the 2006 Eot-Wash 
experiment 

The 2006 Eot-Wash experiment HH searched for de- 
viations from the 1/r 2 force law of gravity. The experi- 
ment used two parallel plates, the detector and attractor, 
which arc separated by a (smallest) distance d = 55/im. 
The plates have holes of different sizes bored into them, 
and the attractor is rotating with an angular velocity w. 
The rotation of the attractor gives rise to a torque on the 
detector, and the setup of the experiment is such that this 
torque vanishes for any force that falls off as 1/r 2 . In be- 
tween the plates there is a d s = 10 /im BcCu-sheet, which 
is for shielding the detector from electrostatic forces. 

The chameleon force between two parallel plates, see 
e.g. Ref. [iH, usually falls off faster than 1/r 2 , imply- 
ing a strong signature on the experiment. However, if 
the matter-coupling is strong enough, the electrostatic 
shield will itself develop a thin-shell. When this hap- 
pens, the effect of this shield is not only to shield elec- 
trostatic forces, but also to shield the chameleon force on 
the detector. This suppression is approximately given by 
a factor exp(— m s d s ), where m s is the mass inside the 
electrostatic shield. Hence the experiment cannot detect 
strongly coupled chameleons. 

The behavior of chameleons in the Eot-Wash experi- 
ment have been explained in Refs. [20|, Hll, HH . We calcu- 
late the Eot-Wash constraints on our models numerically 
based on the prescription presented in Ref. 27 1. 
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FIG. 3: The combined local gravity constraints on the poten- 
tial (|56[1 with p = 1 in the (n, Q) plane. The shaded region 
corresponds to the allowed parameter space. The natural val- 
ues of Q and n of the order of unity are excluded. 



C. Combined local gravity constraints 

1. Potential V(4>) = M 4 exp[p(M '/ '4>) n \ 

Let us first consider the inverse-power law potential 
([56]) with n > 0. From Eq. ([12]) the field value (j> m at the 
minimum of the effective potential V c r satisfies 



M 



n+l 



pn M p i V{(j) m ) 



(70) 



In this model the field is in the regime M <C cf> m <C M p \ 
for the density p m we are interested in. Since V(4> m ) can 
be approximated as V(<fi m ) — M A , it follows that 



M p v 



M \pnM pl M 4 



-l/(n+l) 



(71) 



Using the LLR bound ([41]) with the homogeneous den- 
sity p m ~ 10 _24 g/cm 3 in our galaxy, we obtain the con- 
straint 



n . 10 10-16n < Q/ jJL _ 

The WMAP bound (J69j) gives 



^ M 



' 0M n+1 p£ r 
pn M 4 



i/i, 



(72) 



(73) 



where is the matter density today, with pm' /M 4 sa 

^m O) /^0 O) ~ I/ 3 - Sincc Mpi/M « 10 30 , this condition is 
well satisfied for Q, n, of the order of unity. 



(0) 



For the Eot-Wash experiment, the chameleon torque 
on the detector was found numerically to be larger than 
the experimental bound when Q, n, p are of the order of 
unity. Providing the electrostatic shield with a thin-shell, 
we require that n > 1, Q > 1 or ;t < 1 to satisfy the 
experimental bound. 

In Fig. [3] we plot the region constrained by the bounds 
([72"]). ([75]) . and the Eot-Wash experiment for p, = 1. This 
shows that only the large coupling region with Q ^> 1 can 
be allowed for n of the order of unity. A viable model can 
also be constructed by taking values of p much smaller 
than 1. Note that the WMAP bound ([73]) is satisfied for 
the parameter regime shown in Fig. [3j 



2. Potential V((f>) = M 4 [l - p(l - e^'^y 



Let us proceed to another potential (|58|) with < n < 
1. In the regions of high density where local gravity ex- 
periments are carried out, we have <p <C M v \ and hence 
V{(j>) ~ M 4 [l - p((f>/M p i) n ]. In this regime the effective 
potential V c g has a minimum at 



Q Pm 

pn M 4 



i/(n-i) 



P l ■ 



(74) 



Recall that the LLR bound was already derived in 
Eq. ([rj2"]) , which is the same as the constraint ([72"]) of the 
previous potential. 

Assuming that the chameleon is at the minimum of its 
effective potential in the cosmological background today, 
the WMAP bound (69]) translates into 



Q > 0.05(60n/i) 1 /" 



(75) 



where we have used p^n /M 4 w 1/3 in Eq. flTJ]). However, 
a full numerical simulation of the background evolution 
shows that this is not always the case. For a large range 
of parameters the chameleon has started to lag behind 
the minimum, which again leads to a weaker constraint. 

The Eot-Wash experiment provides the strongest con- 
straints when Q is of the order of unity for the potentials 
([ST]) and ([56]) . This is not the case for the potential (|58| . 
because the electrostatic shield used in the experiment 
develops a thin-shell. 

The mass inside the electrostatic shield is given by 



m 2 s ~ n(l — n)p 



~ M 



Using p s 



^pnM 4 

10g/cm 3 and d s ~ 10 pm we have 



T7T- (76) 



M 2 pl 



n)p 



/'" 



10 _ 



(77) 



Taking Q and p to be of the order of unity, we find 
m s d s > 1 as long as n > 0.2. The suppression of the 
chameleon torque due to the presence of the electrostatic 
shield makes the chameleon invisible in the experiment. 
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FIG. 4: The combined local gravity constraints on the poten- 
tial (|58[) with (i = 0.5 in the (n, Q) plane. In this case the 
allowed parameter space (shaded region in the figure) is de- 
termined by the WMAP constraint and the LLR constraint. 
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FIG. 5: The combined local gravity constraints on the po- 
tential (|58l) with fi = 0.05. The allowed parameter space is 
determined by the LLR bound. The WMAP constraint is 
satisfied for the whole parameter space in the figure. 

In Figs.|4]and[5]we plot the allowed regions constrained 
by the bounds (|62p . (ITS]) , and the Eot-Wash experiments, 
for fi = 0.5 and fi = 0.05, respectively. When p, = 0.5 the 
WMAP constraint gives the tightest bound for Q > 0.2, 
and the parameter space Q < 1 is viable for n > 0.7. If we 
decrease the values of p, down to 0.05, then the WMAP 
bound is well satisfied for the parameter space shown in 
Fig.[U Instead, the LLR experiment provides the tightest 
bound in such cases. When /i = 0.05, the region with 
0-1 ^$ Q ^ 10 and n > 0.6 can be allowed. The coupling 



Q as well as the parameter n are not severely constrained 
for the potential ([58]) . 



VI. LINEAR GROWTH OF MATTER 
PERTURBATIONS 

We now turn our attention to cosmological perturba- 
tions in chameleon cosmology. It is well-known that mat- 
ter perturbations allow to discriminate between dark en- 
ergy models where the gravitational interaction is mod- 
ified on cosmic scales. We first review the general for- 
malism and derive the equation for linear matter pertur- 
bations. We also introduce important quantities like the 
critical scale A c below which modifications of gravity are 
felt and the growth index j(z,k), a powerful discrimina- 
tive quantity for the study of the modified evolution of 
matter perturbations as was explained in the Introduc- 
tion. 

As we have seen in Sec.fVl local gravity constraints im- 
pose very strong boundaries on the potential (|56|) . forc- 
ing its parameters to take unnatural values. Moreover, 
for viable choices of parameters, we have verified that 
the linear perturbations behave in a manner similar to 
the ACDM model, as A c is much smaller than the cosmic 
scales we are interested in. On the other hand we have 
shown that the model parameters of the potential ([55)) 
is not severely constrained. In Sec. IVIBI we will show 
that the potential (|58| gives rise to some very interesting 
observational signatures. 

A. General formalism for cosmological 
perturbations 

We consider scalar metric perturbations a, B, ip, and 
7 around a flat FLRW background. The line-element 
describing such a perturbed Universe is given by |(34[ 

ds 2 = -(1 + 2a)di 2 - 2aB 4 dtdx i 

+a(t) 2 [(1 + 2^)Sij + 2 7)i;i ] dxW • (78) 

We decompose the field <fi into the background and inho- 
mogencous parts: cf>(t, x) = <f>(t) + S(p(t, x). The energy- 
momentum tensors T^y of non-relativistic matter can 
be decomposed as 

n [m) = ~(p* m + 5p* m ) , Tf n) = -p* mVli , (79) 

where v is the peculiar velocity potential of non- 
relativistic matter. In the following, when we express 
background quantities, we drop the tilde for simplicity. 

Let us consider the evolution of matter perturbations, 
S m = Sp* n / ' p* m in the comoving gauge (v = 0). The quan- 
tity S m corresponds to the gauge-invariant quantity in- 
troduced in Refs. [63 - [r36l | when expressed in the comoving 
gauge. In the Fourier space the first-order perturbation 
equations are given by [671 [68j 
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a = -Qty/Mpi , (80) 

5p* m + 3HSp* m - p* m (k - 3Ha) - Q{p*Jcf> + 6p* m j>)/M pl = , (81) 

S<f) + 3H5(f> + (v M + ^jScj) + 2aV t4> - <j)(a - 3Ha + k) + Q{2ap* m + 5p* m )/M pl = , (82) 

k + 2Hk + 3Ha - ^ a - — L- - 4a0 2 + 4^0 - 2V^<5<?(>) = , (83) 



r 



where k is a comoving wavcnumber and k = 3(Ha — ?/>) + where we have used Eq. ()80p . Plugging Eq. (|84p into 
(k 2 /a)(B + a-y). From Eq. dHU) it follows that Eqs. ([82]) and ([83]), we obtain 

K = S m - Q{H + 3HScj))/M p i , (84) 



3H + 2Q 



pi 



v a 2 M p i 



Pin 

Q M k 



8<j>- 



Pm$m 

2M 2 



(1-2Q 2 ) + 



M pl 



_Q_ 

" M pl 

2fc2 

^2" " 



6# 2 - 6# 



<fi $m = , 

20 2 



(85) 



Mpi 



(86) 



where m 2 = V ^ is the mass squared of the chameleon 
field. 

As long as the field 4> evolves slowly ( "adiabatically" ) 
along the instantaneous minima of the effective potential 
V e ft, one can employ the quasi-static approximation on 
sub-horizon scales (k aH) [65], |66], [69| . This corre- 
sponds to the approximation under which the dominant 
terms in Eqs. ([85]) and ([86]) arc those including k 2 /a 2 , 
ui(p , and 5 m , i.e. 

(ml + ^5<j>~--Q-p* m 8 m , (87) 

S m + 2HS m -^^(l-2Q 2 ) 
lM P \ 

+ £ fa + m 

where we have also used the approximation <p/M p \ <C H. 
Combining these equations, it follows that 

5 m + 2H5 m - iTtG eS p* m 8 m ~ , (89) 

where the effective gravitational coupling is given by 

G -^ G ( 1+2<?2 ^w)' (90) 



An analogous modified equation was found in Refs. [ll|, 
UlI j the crucial point being to elucidate the physi- 
cal significance of G c ff. Wc can understand the physi- 
cal content of the modification of gravity by looking at 
the corresponding gravitational potential in real space. 
The gravitational potential (per unit mass) is of the type 
V(r) = ~{G/r) (1 + 2Q 2 e - m * r ) [HI- 

When we solve the full system of perturbations (1861) . 
we can, for some of our models, get a small discrepancy 
compared to ([89]) . This can result in a non- negligible dif- 
ference of up to around 5% in the numerical calculation 
of the growth rate of matter perturbations. This arises 
mainly because the field cj> does not move exactly along 
the minimum of the effective potential but is instead lag- 
ging a little behind it. 

We see that in chameleon models G c g is a scale- 
dependent as well as a time-dependent quantity. Clearly 
the scale-dependent driving force in Eq. (j90p induces in 
turn a scale dependence in the growth of matter pertur- 
bations with two asymptotic regimes, i.e. 

G off = G(l + 2Q 2 ) k/a > m or A « A c , (91) 
= G k/a < m<£ or A > A c , (92) 

where we have introduced the physical wavelength A = 
(2ir/k)a. Wc have in particular Aq = 2n/k today (a = 1). 
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The characteristic (physical) scale A c is defined by 
A c = 2-7t/to0 . 



(93) 



On scales A > A c matter perturbations do not feel the 
fifth force during their growth. On the contrary, on scales 
much smaller than A c they do feel its presence. During 
the matter dominance (fl m ~ 1) the solutions to Eq. ([89]) 
are given by 



S m OC a [V'+M(W)-H/4 

6 m cx a 



A< A c , 
A> A c . 



(94) 
(95) 



Hence, in the regime A <C A c , the growth rate gets larger 
than that in standard General Relativity. 

As mentioned earlier, a powerful way to describe the 
growth of perturbations is by introducing the function 
-f(k, z) defined as follows 



where 



/ = 



d In S m 
din a 



(96) 



(97) 



We remind the definition S m = 5p* n /p* n . The quantity 
7 can be time-dependent and also scale-dependent. It 
is known that a large class of dark energy models in- 
side General Relativity yields a quasi-constant 7 with 
values close to that of the ACDM model, 7 ~ 0.55 
pH l46j ] . Therefore any significant deviation from this 
behavior would give rise to a characteristic signature for 
our chameleon models. Since < Q m (z) < 1, smaller 7 
implies a larger growth rate of matter perturbations. 

As we have seen before, the chameleon mechanism is 
devised so that in high-density environments the mass of 
the scalar field is large relative to its value in low-density 
ones. During the cosmological evolution, the mass of the 
field will follow this behavior, which means that A c will 
move from small to large values. In other words, G c ff 
will evolve from the regime ([92)) to the regime ([91]) . This 
transition is scale-dependent, which is an important fea- 
ture of the growth of matter perturbations in chameleon 
(and f(R)) models. 

We consider the evolution of matter perturbations for 
the wavenumbers 



0.01 /iMpc" 1 < k < 0.2 tiMpc^ 1 



(98) 



where h describes the uncertainty of the Hubble param- 
eter Hq today, i.e. Hq = 100 /ikmsec -1 Mpc -1 . The 
scales (|98|) range from the upper limit of observable scales 
in the linear regime of perturbations to the mildly non- 
linear regime (in which the linear approximation is still 
reasonable) . 

Depending on the value of A c today (denoted as A C; o) 
and on its recent evolution, three possibilities can actu- 
ally arise: (i) The model is hardly distinguishable from 
ACDM; (ii) The model is distinguishable from ACDM 



but shows no dispersion i.e. no scale-dependence. In 
this case low values of 70 will also yield large slopes 
7q = (drf/dz)(z = 0), much larger than in ACDM; and 
finally (iii) The model is distinguishable from ACDM and 
shows some dispersion altogether. These three cases can 
be characterized using the quantity 70 = 7(2 = 0). This 
classification is an alog ous to what was done for some vi- 
able f{R) models [54J and it can be defined as follows: 



• (i) The region in parameter space for which 70 > 

0. 53 for all the scales described by ([9"8|) . In this 
region 7q is small and 7 is nearly constant. 

• (ii) The region where 70 is degenerate, i.e. as- 
sumes the same value for all the scales, with a value 
smaller than 0.5. In this region 7q is large and there 
is a significant variation of 7. 

• (iii) The region where 70 shows some dispersion, 

1. e. a scale-dependence. For low 70, 7g is large and 
we have significant changes of 7. 

Models in the regions (ii) and (iii) can be clearly discrim- 
inated from ACDM. Some examples are shown in Figs. [6] 
and [7] We investigate below in more details the appear- 
ance of these characteristic signatures. 



B. Observational signatures in the growth of 
matter perturbations 

The main question when looking at linear perturba- 
tions in chameleon models is the order of magnitude of 
the scale A c . If the chameleon mass is large such that 
A c is less than the order of the galactic size, matter per- 
turbations on the scales relevant to large scale structures 
do not feel the chameleon's presence. This is actually 
the case for the inverse power exponential potential (|56p 
with model parameters bounded by observational and 
local gravity constraints. Therefore, we will concentrate 
the analysis of the perturbations on the potential (|58jl . 

In the regime 4>/M p \ <C 1 one can employ the approxi- 
mation V((j>) ~ M 4 [l — /u(<^/M p i)™]. In fact this approxi- 
mation is valid for most of the cosmological evolution by 
today. Then we obtain the field mass at the minimum of 
the effective potential V e s{<t>)'- 



2TT 3(2- n ) 

e 2 (!-") 

A c ,o 



N 



(99) 



where A Ci q is the critical length today, given by 



2n 



v.o 



M p i I Q p 



(0) \ 2(l-n) 



vVi(l - n) M 2 \pnM 4 



(100) 



Using the relation p^ 



M , we find 



that A c .o is at most of the order of H n 



If Q 



1. 



n = 0.6, /i = 0.05, and fi^g = 0.72, for example, A c ,o ~ 
OAHq . For the modes deep inside the Hubble radius 
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FIG. 6: The evolution of the growth index 7(2) versus the 
redshift z in the model (|58l) with n = 0.8, /i = 0.5 and 
Q = l/v6 for four different values of k. For these param- 
eters, the model passes local gravity constraints. In this case 
the field is lagging behind the minimum, which means that 
the approximated equation (|89[1 provides the results slightly 
different from this plot (by around 4% in the value of 70). 



today (Ao <C H 1 ) the perturbations are affected by the 
modification of gravity. 

Plugging Eq. ([99]) into Eq. ([90|) . the effective gravita- 
tional coupling can be expressed in the form 



where 



G e ff — G 
A = 



2Q 2 



1 + e-MN-Nc) 



1-n' 
2(1 - n) 
i-n 



In 



A c ,o 



(101) 

(102) 
(103) 



Here A and N c are the parameters describing respectively 
the steepness of transition from the inferior to the supe- 
rior asymptote and the position of the half-amplitude 
value of the function [i.e. G fr = G(l + Q 2 )]. 

The first interesting point to remark is that the steep- 
ness of transition depends exclusively on the value of n, 
with a step function as a limit when n — > 1 and a slower 
transition for n — > 0. Nonetheless, it must be remem- 
bered that this is so in the variable N. When converting 
back to the redshift z, the logarithm scale introduces a 
distortion, which means that the earlier a scale starts its 
transition, the slower this transition will be. In this sense, 
there is a scale-dependence to the steepness of the tran- 
sition, even if it does not appear explicitly in Eq. (|102j) . 

If we want to get a better feel for the epoch of the 
transition, we can define the redshift z c = c~ Nc — 1. From 



FIG. 7: The evolution of 7 versus the redshift z in the model 
(158[) with n — 0.7, /1 = 0.05 and Q — 1 for four different values 
of k. In this case the model passes local gravity constraints, 
too. Contrary to Fig. [6] we have a significant dispersion of 7, 
depending on the wavenumbers k. 



Eq. (fT03"l) it follows that 



A c .o 
^0" 



2(1-") 



(104) 



As expected, the transition redshift is scale-dependent, 
with higher values for the smaller scales. Moreover, 
z c — > for all scales when n — > 1. Since this is also 
the limit at which the transition becomes a step func- 
tion, we have thus an indication that a transition that 
happens very close to the present will necessarily be very 
steep. Another remark is that the transition happens 
in the past (z c > 0) for the scales Ao smaller than A Ci o- 
Thus A Cj0 is a good indication of the scale around which 
there will be a dispersion, as it marks the scale that will 
be exactly in the middle of the transition today. 

If we write z c explicitly in terms of the parameters and 
constants of the model, we obtain 



(/m)' 



1-n M 4 



Q 



o ioy 

M 1 



1. (105) 



This shows that z c gets larger with increasing (i (as the 
deviation from the ACDM model is more significant) 
and/or decreasing Q. Although the transition occurs ear- 
lier for a weaker coupling Q, we need to take into account 
the fact that the growth rate in the regime z < z c is 
smaller for a weaker coupling. 

Let us proceed to the numerical analysis of the growth 
of perturbations. In Fig. [6] we plot the evolution of the 
growth indices 7 for Q = 1/V&, n = 0.8 and /i = 0.5 for a 
number of different wavenumbers within the range spec- 
ified by Eq. (f55|) . Note that these parameters satisfy the 
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=3, 




0.2 0.4 0.6 0.8 



n 



=3. 




0.2 0.4 0.6 0.8 



n 

FIG. 8: The regions (i), (ii) and (iii) in the (n,fj,) space for 
Q — l/i/E (top) and for Q = 1 (bottom). In the region (i) all 
the modes have 70 > 0.53. In region (ii) 70 is degenerate with 
the value smaller than 0.5, and finally (iii) shows the regime 
where 70 is dispersed. It is clear that, for both choices of Q, 
there are viable choices of parameters in which the deviation 
from the ACDM model is present. 



local gravity constraints discussed in Sec. El see Fig. [4] 
At the present epoch these modes are in the regime (ii), 
with very similar growth indices today (70 — 0.46). As 
estimated by Eq. f| 105f> the transition rcdshift z c is larger 
than the order of 1 for k > 0.05 /iMpc" 1 , e.g. z c = 2.3 
for k = 0.1 /iMpc -1 . The degenerate behavior similar to 
that shown in Fig. [6] has been also found in some f(R) 
models 0, [54| . Numerically we have verified that values 
of 70 in this case are slightly higher than those expected 
in the asymptotic regime (A <C A c ). This discrepancy 
comes from the fact that the chameleon field is lagging 
behind the minimum of the potential for this choice of 



parameters. As a result, we need to solve the full per- 
turbation equations (|85|) and (|86| instead of the approx- 
imated equation ([5§]l . Our numerical results show that 
there can be a discrepancy of up to a few percent in the 
growth rate calculated with the two different methods. 

Another choice of parameters, which are compatible 
with local gravity constraints, is made in Fig. where 
we see the evolution of the growth indices 7 for Q = 1, 
n = 0.7 and \i = 0.05. The behavior of the growth indices 
is very different from the one shown in Fig.[6l Clearly this 
corresponds to the regime (iii), in which 70 is dispersed 
with respect to the wavenumbers k. The reason for the 
dispersion is that the transition to the regime A <C A c 
occurs on lower rcdshifts than in the case shown in Fig. [6] 

In Fig. [8] we illustrate three different regimes in the 
(n, (J.) plane for the couplings Q = l/y/6 and Q = 1. For 
the model parameters close to (n, fi) = (1,0), the pertur- 
bations behave similarly to those in the ACDM model 
[i.e., in the region (i)]. Figure [S] shows that the limits 
imposed by the constraints derived in Sec. El although 
strong, allow the large parameter space for the existence 
of an enhanced growth of matter perturbations and for 
the presence of the dispersion of 70. Finally we note 
that the large variation of 7 in the regime z < 1 seen 
in Figs. [6] and [7] will also enable us to distinguish the 
chameleon models from the ACDM model. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have studied observational signatures 
of a chameleon scalar field coupled to non-relativistic 
matter. If the chameleon field is responsible for the 
late-time cosmic acceleration, the field potentials need 
to be consistent with the small energy scale of dark en- 
ergy as well as local gravity constraints. Wc showed 
that the inverse power potential cannot satisfy both cos- 
mological and local gravity constraints. In general, we 
require that the chameleon potentials are of the form 
V((f>) = M 4 [l + /(</>)], where the function /(</>) is smaller 
than 1 today and M is a mass that corresponds to the 
dark energy scale (M ~ 10~ 12 GcV). 

The potential V(<j>) = M 4 exp[^(M/^) n ] is one of those 
viable candidates. However we showed that the allowed 
model parameter space is tightly constrained by the 2006 
Eot-Wash experiment. As we see in Fig. [31 the natural 
parameters with n and Q of the order of unity are ex- 
cluded for \i = 1. Unless we choose unnatural values of 
/1 smaller than 10 -5 , this potential is incompatible with 
local gravity constraints for {n, Q} = 0(1). 

On the other hand, the novel chameleon potential 
V((j>) = Af 4 [l - ^(1 - e-^p 1 )™], which has the asymp- 
totic form V((f>) ~ M 4 [l — /j,(<j> / M p \) n ] in the regime 
4> <C Mpi, can be consistent with a number of local grav- 
ity experiments as well as cosmological constraints. In 
fact this case covers the viable potentials of f(R) dark 
energy models in the Einstein frame. The allowed param- 
eter regions in the (n, Q) plane are illustrated in Figs. 0] 
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and [5] for /j, = 0.5 and fi = 0.05. This potential is viable 
for natural model parameters and for the coupling Q of 
the order of unity. 

In order to distinguish the chameleon models from the 
ACDM model at the background level, we discussed the 
evolution of the statefmders (r, s) defined in Eq. (f6"3")l . 
Unlike the ACDM model in which r and s arc con- 
stant (r = 1, s = 0) the statefmders exhibit a pecu- 
liar evolution, as plotted in Fig. [2] For the potential 
V{4>) = M 4 [l - fji(l - e -*/ M Pi)«] we found that r < 1 
and s > around the present epoch, but the solutions 
approach the de Sitter point (r, s) = (1, 0) in future. The 
upcoming observations of SN la may discriminate such 
an evolution from other dark energy models. 

We have also studied the growth of matter perturba- 
tions for the chameleon potential V{<j>) — M 4 [l — fi(l — 
e -<P/M P iyi^ The presence of a fifth force between the field 
and non-relativistic matter (dark matter/baryons) mod- 
ifies the equation of matter perturbations, provided that 
the field mass is smaller than the physical wavenum- 
ber k/a, or A < A c . Cosmologically the field is heavy 
in the past (i.e. for large density), but the mass de- 
creases by today (typically of the order of Hq) in order to 
realize the late-time cosmic acceleration. Then the tran- 
sition from the regime k/a < m$ (A > A c ) to the regime 
k/a > m§ (A < A c ) can occur at the redshift z c given in 
Eq. f|105|) . For the perturbations on smaller scales (i.e. 
larger k) the critical redshift z c tends to be larger. 

For the model parameters and the coupling Q bounded 
by a number of experimental and cosmological con- 
straints, we have studied the evolution of the growth in- 
dex 7 of matter perturbations. Apart from the "General 
relativistic regime" in which two parameters n and (i of 
the potential ([58]) are close to (n, /i) = (1,0), we found 
that the values of 7 today exhibit either dispersion with 
respect to the wavenumbers k (region (iii) in Fig. [8]) or 
no dispersion, however with 70 smaller than 0.5 (region 
(ii) in Fig. [8|) . Both cases can be distinguished from the 



ACDM model (where 7 ~ 0.55). Moreover, as seen in 
Figs. [6] and [TJ the variation of 7 on low rcdshifts is sig- 
nificant. 

From observations of galaxy clustering we have not yet 
obtained the accurate evolution of 7. This is linked to 
the fact that all probes of clustering are plagued by a 
bias problem. However upcoming galaxy surveys may pin 
down the matter power spectrum to exquisite accuracy, 
together with a better understanding of bias. In order 
to confirm or rule out models like ours, one must also 
address the observability of 7 both function of z 

and k. We hope that future observations will provide an 
exciting possibility to detect the fifth force induced by 
the chameleon scalar field. 
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